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The behaviour of a power transformer is complex and difficult to predict during tran- 
sient conditions or during operation at frequencies below or above its nominal fre- 
quency, a phenomenon common in renewable energy plants due to harmonic distor- 
tion. Furthermore, the accuracy of a power system simulation depends on the models 
of critical subsystems such as the power transformers. This paper presents the use of 
a unique excitation waveform comprising of pseudo-random current impulses to ac- 
curately identify the wideband characteristics of a power transformer. By injecting 
the excitation waveform to the relevant transformer terminals, frequency responses are 
determined by cross-correlation of the perturbation signal, and the measured response. 
Compared to the traditional transformer identification methods, the pseudo-random 
current impulses offer a wideband excitation with a higher degree of controllability 
such that its spectral energy can be focused in the frequency band of interest. The pro- 
posed method was investigated on a 16 kVA, 22 kV/240 V single-phase transformer. 
The obtained wideband frequency responses provide useful information in harmonic 
penetration and over-voltage studies and are also used to estimate, with a high degree 
of accuracy, the lumped parameters of the equivalent transformer broadband circuit 
model. 
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1. INTRODUCTION 


In recent years, there has been a rapid increase in renewable energy sources such as photovoltaic and 
wind generators [1]-[3]. These energy sources have presented new challenges such as the risk of high harmonic 
voltage distortion due to their use of power electronic converters [1], [4], [5]. These power system harmonics 
have increased in both amplitude and frequency bandwidth. Step-up transformers that are used to connect these 
renewable plants to the grid are adversely affected by the harmonics from these renewable plants [6]. 


Power transformers do not have a flat frequency response, which makes them unsuitable for use in 
applications with high harmonic components and high frequency transient signals. Transformer impedances 
are functions of the input signal frequencies and voltage, as well as the loading conditions [7], [8]. Transformer 
parameters such as leakage inductances, winding capacitances, and the non-linear core components interact 
to cause complex series and/or parallel resonant phenomena in the transformer frequency characteristics [9]. 
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When a minimally damped resonant point in the power transformer’s wideband transformation ratio charac- 
teristics is excited by a harmonic voltage or a transient event, excessive voltage and/or current amplification 
occurs. Over-voltages and over-currents may cause protection devices to activate, or cause damage to the power 
transformer and other equipment connected to the network, resulting to unnecessary downtime. 

In this context, the response of transformers to input signals with multiple frequencies needs to be 
modeled and understood. Therefore, this gives rise to renewed interest in performing broad frequency band 
measurements, especially on power transformers that are connected to wind and photovoltaic (PV) energy 
sources. A common approach in transformer identification is the use of non-parametric methods that compute 
the transformer frequency responses. Frequency response analysis (FRA) has been widely used for condition 
monitoring of transformers [10]-[15]. Frequency response measurements have also been applied in transformer 
parameter estimation [10], [16]-[18]. Although power transformer frequency responses have been widely used 
for condition monitoring and parameter estimation [10]-[12], [18] their use in power quality studies, especially 
in cases where the transformers are used in renewable energy sources, has not been given much attention. 

In practice, frequency response measurements are conducted by applying an excitation signal that can 
perturb the target system at the frequencies of interest. The system input and output signals are measured over 
an interval for further processing to obtain the relevant frequency responses. The excitation signals proposed 
in literature for transformer FRA include the swept sine and the impulse [11], [19], [20]. The sine sweep 
involves a variable-frequency sinusoidal signal which is simple to implement using low-voltage sinusoidal 
signal generators. In practice, however, it is difficult to design and implement a sinusoidal source to perturb a 
large power transformer over the frequency range of interest, especially at low and high ends of the range. This 
is due to reduced magnetizing reactance and winding capacitive reactance at these frequencies respectively. 
Consequently, very low signal levels result in measurement errors due to a low signal to noise ratio (SNR). 
The impulse method involves perturbing the transformer with a classical voltage impulse Which has a low rise 
time and a higher fall time. Although this method is more accurate, the spectral energy of the classical impulse 
waveform drops rapidly with frequency which would result to limited accuracy in some frequencies. 

Controlled signals have been proposed for use in power system equipment identification [21]-[23], and 
have advantages in that their spectral characteristics can be controlled and are repeatable. Therefore, accurate 
measurements are possible at the expense of complex instrumentation [11]. Despite the advantages presented 
by controlled perturbation signal method, there has been limited application of these signals for transformer 
FRA. In the case of a pseudorandom binary sequence (PRBS), the clock frequency and sequence length can 
be controlled to induce persistent excitation for the dynamic modes associated with the device under test. The 
PRBS has a flat frequency spectrum at frequencies below the -3 dB point which allows for perturbation with 
uniform energy [24]. The PRBS is however not suitable as a perturbation signal for transformers due to a 
possibility of injecting low frequency components that would cause core saturation. This particularly occurs 
when a long PRBS is used [25]. The main contribution of this work, therefore, includes the following: 

— The use of a wideband controlled perturbation signal to obtain accurate transformer frequency responses. 
— The use of transformer frequency responses for power quality studies and parameter estimation. 

In this paper the use of a unique signal consisting of bipolar pseudo-random current impulses as a 
broadband excitation signal to determine the frequency responses of a transformer in the harmonic frequency 
range is proposed. The obtained responses are then used to study the effects of harmonics especially on the 
voltage transfer function, to characterise the power transformer over a wide frequency band as well as to 
obtain transformer wideband parameters. The sequence combines the advantages of the PRBS and the classical 
double exponential type impulse. An advantage of using this signal is that its power spectrum can be controlled, 
for instance, in reducing the energy in low frequencies to avoid transformer saturation during measurements. 
Properties of the signal such as the clock frequency and time constants can be used to focus the spectral energy 
to the frequency range of interest and improve the SNR. The bipolar nature of the sequence ensures that the 
transformer is not driven towards a biased offset point from the normal operating point whereas its wideband 
nature allows tests to be conducted within short time periods compared to conventional methods such as the 
sine sweep. 

In section 2, a methodology for transformer frequency response measurements is discussed. Experi- 
mental results and parameter estimation of a high frequency power transformer model using the experimental 
frequency responses is discussed in section 3. The conclusion is presented in section 4. 
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2. RESEARCH METHOD 
2.1. Transformer testing using a sequence of bipolar pseudo-random current impulses 

The idea of using a bipolar pseudo-random impulse sequence as a perturbation signal in the energy 
systems and power equipment was first proposed in [26]. The bipolar pseudo-random impulse sequence is a 
wideband signal consisting of a series of chopped positive and negative impulses. In this work, a sequence of 
bipolar pseudo-random current impulses was used as a perturbation signal for transformer identification and 
modelling. Figure | illustrates the bipolar pseudo-random impulse sequence. Bipolar pseudo-random current 
impulses can be generated using an efficient circuit topology consisting of a DC source, a full-bridge whose 
switches are controlled by a PRBS gate signal and a series RLC network as shown in Figure 2. 


PRIS (t) 
T 
L 


fi 
Time 


Figure 1. Bipolar pseudo-random impulse sequence 


PRBS Generator 


Figure 2. Circuit diagram for the bipolar Pseudo-random current impulses generator 


The PRBS is a random, deterministic and periodic signal that can be produced using shift registers 
driven by a clock sequence. The occurrence and duration of the binary states is random. The data length for 
one PRBS period is defined as N = 2” — 1 bits, where n denotes the number of shift registers. The series 
RLC network provides a means of optimizing the time-constants associated with the rise-and fall-times of the 
bipolar Pseudo-random current impulses. The circuit of Figure 2 can be viewed as an application of a step 
input to a series RLC circuit, where the step is controlled by the PRBS signal. The sequence of bipolar current 
impulses generated by the circuit in Figure 2 over one period can be expressed as: 


eet len (Yo — Vac + 2£) — e*/2 (Vo — Vde + 1t)] T1T2 


i(t) = = Ti T2 


i=0,2,4... L(T2 = T1) 
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7 a. L(t2 — 71) * [e-t -u(t = tis1)]. 
(1) 


where u(t —t;) — u(t —t;41) depict random switching intervals in the full bridge circuit, 7, and T2 are the time 
constants associated with the series RLC network, Jp and Vo denote the initial current through the inductor and 
the initial voltage across the capacitor respectively. 

In (1) indicates that the excitation current depends on the DC source, the RLC network time constants 
and the PRBS switching intervals which in turn are dependent on the PRBS clock frequency. In practice, 
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the shape of the perturbation signal is also affected by the initial energy in the inductor and capacitor in each 
switching cycle. The properties of the current impulses can easily be controlled by varying excitation source 
parameters namely the input DC source, PRBS clock frequency fe, and the impulses time constants. This 
high degree of control offers capability to improve the signal’s SNR, guarantee repeatable experiments and 
allow the user to concentrate the signal’s energy in the required frequency band. The broadband nature of the 
sequence ensures that the transformer identification procedure is conducted fast compared to using conventional 
excitation signals such as stepped sine. This is an advantage especially in in situ power system measurements 
where parameters such as harmonics are time-dependent. The bipolar pseudo-random impulse sequence can, 
therefore, be adjusted to achieve persistent excitation for a broad range of frequencies and applications. This 
is an important aspect especially in identification of systems involving electromagnetic components such as 
transformers where superior performance is achieved when excitation energy is directed in the higher frequency 
band and the low frequency perturbation is restricted. 


2.2. Transformer test arrangement 
2.2.1. Offline frequency response measurements 

In offline frequency response measurements, the transformer is disconnected from the normal AC 
supply voltage. A perturbation source is then connected to the transformer terminals thus providing excitation. 
Figure 3 shows a test arrangement that can be used to obtain the frequency response of the secondary open- 
circuit input impedance. This is the impedance seen from the secondary terminals of transformer T, with 
primary terminals left open. A perturbation source is connected to the secondary terminals as shown in Figure 
3. The waveforms of the applied perturbation signal, i,(t), and the perturbed voltage across the secondary 
winding, v,(t), are then recorded over a time interval. 


vp(t) | E vt) 


(0; 


Figure 3. Test arrangement for secondary open-circuit input impedance frequency response 


Regarding the transformer as a linear time-invariant system and ignoring measurement and quantiza- 
tion noise, the sampled voltage is given by the discrete convolution of the impulse response with the perturba- 
tion signal which can be expressed as. 


vsin] = X` z[k]ip[n — k] (2) 
k 


al 


where v, (n) is the discrete-time voltage signal, ip(k) is the discrete-time perturbation signal and z(k) is the 
sampled system impulse response. 

Convolution of time domain signals implies multiplication of the signals in frequency-domain and 
therefore, the secondary open-circuit input impedance frequency response Z, (w) can be obtained by applying 
the Fourier transform to the measured voltage and current signals. Zs (w) can then be expressed as. 


zlo) = Piet} 


= Fli,()} = 


where F denotes the Fourier transform. 

The primary open-circuit input impedance frequency response Zp „„e„ (w) can be obtained from the 
test arrangement shown in Figure 4. The excitation source is connected to the primary winding terminals and 
the perturbation current, ip(t), and the perturbed voltage across the primary winding, vp(t), are then recorded 
over a time interval. 
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Figure 4. Test arrangement for primary open-circuit input impedance frequency response 


The frequency response of the primary open-circuit input impedance can be obtained from the expres- 


F{vp(t 
AR a @) 
F{ip(t)} 
The frequency response of the transformation ratio TpR(w) can be obtained using the topology shown in [4] 
by recording the perturbed voltages across both primary and secondary windings. The transformation ratio 


frequency response can be expressed as. 


sion 


_ Fup (t)} 
Fiva(t)} 


By short-circuiting the secondary terminals of the test arrangement shown in Figure 4 the primary short-circuit 
input impedance frequency response can be obtained using the expression in (4). 


Tr(w) (5) 


2.2.2. In situ frequency response measurements 

In in situ measurements, the transformer stays connected to the conventional supply voltage during 
perturbation. The transformer does not need to be disconnected, which reduces down time. Furthermore, 
during in situ measurements, the effects of the normal transformer operation are accounted for [26]. The 
dynamic range of the excitation signal however has to be kept low enough to ensure that the transformer is not 
perturbed too far from its optimal operating point. 

Figure 5 shows a connection scheme that can be used to characterize, in situ, the secondary open- 
circuit input impedance Zs. The transformer is connected to its normal ac voltage supply Vra which has 
an internal impedance Zr. The perturbation source is then connected across the voltage source, thereby 
injecting the excitation current into the transformer terminals. The perturbation current, ip(t), causes a voltage 
perturbation to be superimposed on the supply voltage Vra due to a voltage drop across Zp. 


f Z a 
is(t) Th iin (t) 
o Tr <4 C4 <4 


E p ip(t) 


Up(t) 3 f Us(t) ©) vra (t) 


Perturbation source 


Figure 5. Test connection for in situ transformer perturbation 


The perturbed voltage vs(t) and current i,(t) into the transformer are then measured over a time 
interval. The secondary open-circuit input impedance frequency response Z,(w) can be obtained from the 
expression, 


Pais ann (6) 


A similar arrangement can be used to obtain the voltage transformation ratio whereas the secondary 


short-circuit input impedance measurement can be done with the primary-side terminals shorted. In in situ 
measurements, a high SNR can be achieved with lower excitation compared to the offline measurements. 
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3. RESULTS AND ANALYSIS 
3.1. Frequency response analysis 

Verification of the proposed perturbation method is done by perturbing, offline, a 16 KVA, 22 kV/240 
V single-phase distribution transformer using a pseudo-random current impulse described in section 2.1. The 
perturbation signal is generated by controlling the H-bridge switches using a PRBS gate signal with f, set to 
15 kHz thereby switching a series RLC network with values R = 50, L = 100wH,C = 2uF. A DC voltage 
of 60 V is applied to the H-bridge input. 

The excitation current generated is then applied to the transformer terminals using the test arrange- 
ments discussed in section 2.2 to obtain the various frequency response estimates. Each measurement is con- 
ducted for a duration of 1s and the required voltage and current data recorded. Voltage and current waveforms 
are recorded using the GWinstek GDP-025 differential voltage probe and the A622 Tektronix current probe 
respectively. The National Instruments CompactDAQ chassis and the NI 9223 module is used for data acqui- 
sition. Figure 6 shows the measured waveforms for the bipolar Pseudo-random perturbation current and the 
PRBS gate signal obtained for the circuit shown in Figure 2. 


Figure 6. Experimental waveforms for the bipolar Pseudo-random perturbation current (top) and PRBS gate 
signal (bottom) 


To eliminate stochastic measurement noise and reduce leakage and aliasing problems, the acquired 
data is processed in Matlab using the Welch estimator where the measured data is divided into overlapped 
and windowed segments. The transfer function of each segment is computed and then averaged to generate 
a frequency response. The Hann window was selected as it offered good frequency resolution and an overlap 
percentage of 50% is used. The low voltage (LV) open circuit input impedance is measured by using the test 
arrangement shown in Figure 3. 

In this case, the bipolar pseudo-random current impulses generator is connected to the LV side termi- 
nals of the transformer while the high voltage (HV) winding is left open. The perturbed LV Voltage and the 
perturbation current are recorded and used for the impedance estimation. The obtained experimental LV open 
circuit input impedance magnitude and phase frequency responses are shown in Figure 7. 

The impedance frequency response in Figure 7 indicates various resonance peaks. The first resonance 
peak indicates a parallel resonance at approximately 100 Hz, a series resonance at about 6 kHz and a parallel 
and series resonance pair between 20 kHz and 30 kHz. A highly damped parallel and series resonance pair is 
also noticeable between 40 kHz and 50 kHz. 

Experiments to determine the transformer input response from the HV side are conducted using the 
test arrangement shown in Figure 4. To determine the HV short-circuit input impedance, the perturbation 
current is applied to the HV terminals with the LV terminals shorted. The perturbed HV Voltage and the 
injected current are recorded and used in the estimation of the impedance response. In this case, a large current 
is drawn by the winding resistances and leakage inductances due to the shorted LV winding. Figure 8 shows the 
obtained experimental short-circuit input impedance magnitude and phase frequency responses of the HV side 
of the transformer. A parallel resonance at about 2.5 kHz and a series resonance at about 3 kHz are evident. 
At low frequencies, up to about 2 kHz, the impedance is predominantly inductive which can be attributed to 
the leakage inductance of the transformer. However, at higher frequencies beyond 4 kHz, the impedance turns 
capacitive which can be attributed to the large HV winding capacitor. 

The transformation ratio is obtained by applying the pseudo-random current impulse sequence to 
the HV side with the LV winding open. Voltage waveforms vyv (t) and vzy (t) are recorded and used to 
estimate the transformation ratio frequency response. Figure 9 shows the obtained experimental HV to LV 
transformation ratio magnitude and phase frequency responses. 
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Figure 7. Experimental LV open-circuit input impedance magnitude (top) and phase (bottom) frequency 
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Figure 8. Experimental HV short-circuit input impedance magnitude (top) and phase (bottom) frequency 
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Figure 9. Experimental HV-LV voltage transformation ratio magnitude (top) and phase (bottom) frequency 
responses 


From the magnitude response of the transformation ratio, it can be seen that at low frequencies up to 
approximately 10 kHz, the ratio Vav (f)/Viv (f) is flat at a value of about 91 as expected. As the frequency 
increases beyond 10 kHz, resonant peak pairs are observed to appear at frequencies beyond 20 kHz. These 
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peaks in the voltage transformation ratio are attributed to the resonance modes caused by the interaction of 
various transformer elements. It is seen that the voltage transfer at 20 kHz is about twice the voltage transfer at 
the fundamental frequency of 50 Hz. The presence of resonance characteristics in the wideband transformation 
ratio such as the one shown in Figure 9 presents challenges to the use of classical distribution transformers 
in applications that are prone to harmonic distortion and high frequency transient signals such as renewable 
energy plants. 


3.2. Transformer equivalent model and parameter estimation 

A transformer model with parameters representing its physical behaviour is important for designers 
[27]. Such a model can also be used for power system simulation. The tested transformer was modelled using 
the circuit presented in Figure 10 which shows the lumped-parameter equivalent circuit model. The circuit has 
been simplified by referring the secondary circuit parameters to the primary side of the transformer. 


is(t) ZTh iin (t) 
o T, < m ` 


t a A ip(t) 


wt) 3E vW O vrli) 
_ _ [] Perturbation source 


Figure 10. Lumped parameter circuit transformer model 


Cp» Cs and Cps denote the transformer primary winding, secondary winding and inter-winding capac- 
itances respectively. R, and Re denote the transformer winding resistance and core losses resistor respectively 
while Ls and Lm denote the leakage inductance and magnetizing inductance respectively. The parameter 
estimation process is summarised in the flowchart shown in Figure 11. 


Populate transformer model 
with initial parameters, 
0 initial 


Estimate the model Obtain transformer 
impedance frequency experimental frequency 
response response 


Zsm(©) 


Objective function 
calculation, €(0) 


Optimization algorithm 


Figure 11. Flowchart of the transformer parameter estimation procedure used in this study 
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The current signal, ip(t), that was measured in the LV open circuit input impedance experiment is 
transformed to the intermediate voltage level using the ratio Te: The transformed current signal is applied 
to a Simulink model of the circuit in Figure 10, whose parameters have been initialised with values based 
on a priori knowledge of the transformer under investigation. The secondary voltage signal at intermediate 
voltage level is then simulated. The transformed current and the simulated voltage signals are used to estimate 
the frequency response of the secondary open-circuit input impedance, Zse(w). The measured and estimated 


responses, Zsm(w) and Zse(w), were compared to yield an objective function. 


The objective function used in this estimation can be stated as. 


M 
e(0) = 5 |Zsm(we) — Zselwr)|?, (7) 
k=1 


where M indicates the number of discrete frequency points in the compared frequency responses and e(0) is 
the objective function and 6 is the vector of parameters to be estimated. The considered equivalent circuit of a 
single-phase power transformer leads to a parameter vector that can be represented as. 


0 = [C,, Rs, Ls, Cs, Re, Lm, Cys (3) 


The optimization procedure was used to minimize the objective function by repetitively changing the trans- 
former model parameters, such that the difference between Zsm(w) and Zse(w) was at its minimum value. The 
jmincon optimization procedure was used as the optimization algorithm to reduce the objective function in this 
investigation. Table 1 presented the estimated parameter values obtained for the test transformer. 


Table 1. Estimated parameter values for the transformer model in Figure 10 


Parameter Value Parameter Value Parameter Value 
Cp 2 nF Cos 0.1375 nF Cs 31.75 pF 
Rs 1k Re 1M 
Ls 0.32 H Lm 1350 H 


The measured and estimated frequency responses of the LV open-circuit input impedance are shown 
in Figure 12. It can be seen that the considered transformer model provides a reasonable good prediction of the 
frequency response up to 20 kHz, despite its relative simplicity. The location of the major resonant frequencies 
of the response, at 100 kHz and 6.2 kHz, and the amount of damping associated with these resonances is well 
predicted. However, the order of the model is too low to accurately predict the minor resonances occurring at 
higher frequencies. The transformer model can be improved by taking into account the non-linear phenomena 
such as magnetic saturation and frequency-dependence of the winding resistances due to skin effect. The 
inductance of the measurement cables between the perturbation source and the transformer terminals, which 
resonates with the transformer stray capacitors, also need to be taken into account in the model. 


Frequencies at which the resonance occurs can be confirmed using (9) 


1 


m TiC 


(9) 


where f, is the resonant frequency, L and C’ denote the inductance and capacitance respectively. 


For instance, the first parallel resonance in Figure 12 occurs when Lm resonates with Cp while the 
first series resonance can be attributed to L, resonating with Cp. To validate the obtained parameters, the trans- 
former model shown in Figure 10 was populated with the parameters in Table 1 and simulated by perturbing it 
using a sine sweep. The model was used to estimate the frequency response of the transformation ratio. Figure 
13 shows that the simulated results match with the measured results in Figure 9 indicating resonance points 
between 20-30 kHz. 


Characterizing power transformer frequency responses using bipolar ... (Fredrick Mwaniki) 


2432 o ISSN: 2088-8694 


Magnitude [Ohms] 


Frequency [Hz] 


Figure 12. Measured and estimated frequency responses of the LV open-circuit input impedance; 1) 
solid-measured, 2) dashed-estimated 
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Figure 13. Transformation ratio frequency response simulated using the estimated parameters 


4. CONCLUSION 

A novel methodology that can be used to study the behaviour of power transformers not only at the fun- 
damental frequency but also at harmonic frequencies is presented. Accurate frequency responses are obtained 
by perturbing a 16 kVA, 22 kV/240 V power transformer using a sequence of bipolar pseudo-random current 
impulses. The time-and frequency-domain properties of this wideband excitation signal are controllable which 
makes it particularly suitable in transformer identification where excitation energy should be concentrated to- 
wards the higher frequency band with limited low frequency components. This reduces the saturation effects of 
the electromagnetic system and is a major advantage of the proposed perturbation signal for transformer iden- 
tification, compared to other approaches such as those utilising sine sweep, impulse and PRBS. The proposed 
method can be used to study the impact of harmonics especially on transformers connected to renewable energy 
sources. Furthermore, the measured frequency response results can be used to estimate the parameter values 
of a wideband lumped parameter equivalent circuit model of the transformer. As future work, the accuracy of 
the estimated wideband transformer parameters can be improved by considering non-linear phenomena such as 
core saturation and skin effect. The non-linearity will result in additional model parameters and hence a higher 
order transformer model. 
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